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Abstract 


We study the stability in the inverse problem of determining the time dependent zeroth-order 
coefficient q{t, x) arising in the wave equation, from boundary observations. We derive, in dimension 
n >2 , a log-type stability estimate in the determination of q from the Dirichlet-to-Neumann map, in 
a subset of our domain assuming that it is known outside this subset. Moreover, we prove that we can 
extend this result to the determination of q in a larger region, and then in the whole domain provided 
that we have much more data. 
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1 Introduction 

1.1 Statement of the problem 

Let be a bounded domain of M”, n > 2, with C°° boundary L = dQ. Given T > 2Diam(r2), we 
introduce the following initial boundary value problem for the wave equation 



( 1 . 1 ) 


where / G uq G ui G L‘^{Q) and the potential q G C^{Q) is assumed to be real valued. 

It is well-known (see ifTSl . lISl ) that if the compatibility condition is satisfied, then ( |1.1| ) is well-posed. 
Therefore we can introduce the following operator 



usually called the Dirichlet-to-Neumann map. Here zz(x) denotes the unit outward normal to L at x and 
dpU stands for Vu.u. 
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In the present paper, we will first study the inverse problem of recovering the time dependent po¬ 
tential q from the Dirichlet-to-Neumann map Aq associated to the problem ( |1.1| with (uq, ui) = (0,0). 
This inverse problem is to know whether the knowledge of Ag, can uniquely determine the electric time 
dependent potential g. 

Physically, it consists in determining physical properties such as the time evolving density of an in¬ 
homogeneous medium by probing it with disturbances generated on the boundary. And the goal is to 
recover q which describes the property of the medium. We assume that the medium is quiet initially and 
the Dirichlet data / is a disturbance used to probe it. 

The problem of recovering coefficients for hyperbolic equations from boundary measurements was 
treated by many authors. In ifTSl Rakesh and Symes proved a uniqueness result in recovering the time 
independent potential q(x) in the wave equation. In ||9l Isakov treated the inverse problem of recovering 
a zeroth order coefficient and a damping coefficient. These results are concerned in the case where the 
Dirichlet-to-Neumann map is considered in the whole boundary. A key ingredient in the existing results, 
is the construction of complex geometric optics solutions concentrating near lines with any direction 
uj G and the relationship between the hyperbolic Dirichlet-to-Neumann map and the X-ray trans¬ 
form plays a crucial role. The uniqueness in the determination of time independent potential appearing 
in the wave equation by a local observations was proved by Eskin f7]. 

The uniqueness by local measurements is solved well. However, the stability by a local Dirichlet-to- 
Neumann map is not discussed comprehensively. For it, one can see Bellassoued, Chouli and Yamamoto 
[3]] where a log-type stability estimate was proved in the case where the Neumann data are observed in 
an arbitrary subdomain of the boundary, Isakov and Sun ifTTl where a local Dirichlet-to-Neumann map 
yields an Holder stability result in determining a coefficient in a subdomain. The case where the Neu¬ 
mann data are observed in the whole boundary, a stability of Holder type was established in Cipolatti and 
Lopez lO, Sun |[23ll . and in Riemannian case in M. Bellassoued and D. Dos Santos Ferreira [4], Stefanov 
and Uhlmann ll2T1l . 

All the above mentioned results are concerned only with time-independent coefficients. Many au¬ 
thors considered the problem of determining time-dependent coefficients for hyperbolic equations. In 
{22} . Stefanov proved that the time dependent potential q{t, x) arising in the wave equation is uniquely 
determined from the knowledge of scattering data. In lfT9l . Ramm and Sjdstrand treated the problem of 
determining the time-dependent potential q{t, x) from Dirichlet-to-Neumann map, on the infinite time- 
space cylindrical domain x D, and they proved a uniqueness result under suitable assumptions. In 
|[20ll . R. Salazar, extended the results in |fT9ll to more general coefficients and proved a result of stability 
for compactly supported coefficients provided T is sufficiently large. 

The inverse problem of determining the time-dependent coefficient q{t, x) from the Dirichlet-to- 
Neumann map Aq, was treated by Ramm and Rakesh |[T6l . they assumed without loss of generality that 
D is a ball and they proved a uniqueness result only in a subset made of lines making 45° with the t-axis 
and meeting the planes f = 0 and t = T outside Q, provided that it’s known outside this subset. It’s 
clear that with zero initial data one can not hope to recover q{t, x) over the whole domain Q, even from 
the knowledge of the full boundary operator Aq. This is due to the domain of dependence associated 
to the hyperbolic problem ( |1.1| ) (see [8]). However, in Isakov lITOl . the ideas from ifTTl - llTSl are used to 
prove a uniqueness result in determining q{t, x) over the whole domain Q, but he needed much more 
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information. Indeed his data was the response of the medium for all possible initial data. 

In this paper, we will prove a log-type stability estimate which establishes that the time dependent 
potential q{t, x) depends stably on the Dirichlet-to Neumann map Aq in a subset of our domain, provided 
that it is known outside this subset. After that we prove that we can extend this result to the determination 
of (7 in a larger region if we further know the measures {u{T ,.), dtu{T ,.)), where u is the solution of the 
initial boundary value problem o> with (uq, ui) = (0, 0). Moreover, we will prove that if our data was 
the response of the medium for all possible initial data, then we have a log-type stability estimate for this 
problem over the whole domain Q. 

Inspired by the work of M. Bellassoued and D. Dos Santos Ferreira |4i, Alden Waters Il24ll succeeded 
in proving a type of an Holder stability estimate for the inverse problem of recovering the X-ray trans¬ 
form of the time-dependent potential g, appearing in the wave equation, from the dynamical Dicrichlet-to 
Neumann map in Riemannian case. A key ingredient in this result is the construction of Gaussian beam 
solutions. In the case n > 3, the inverse problem associated to the system ( |1.1| ) with the initial condition 
uq = 0, was treated recently by Y. Kian |[T2l, indeed, inspired by Bellassoued-Jellali-Yamammoto Gl-lfTl 
and using suitable complex geometric optics solutions and Carleman estimate, he proved a log-log type 
stability estimate in determining the time dependent coefficient q(t,x), from the knowledge of partial 
Dirichlet-to-Neumann measurement and the measure u{T ,.). 

Before stating our main results, we recall the following Lemma on the unique existence of a solution 
to the problem ([U). The proof is given in lfT3l (see also f5]). 

Lemma 1.1 LetT > d be given. Suppose that uq G ui G L^(D), and f G Assume, 

in addition, that /(O,.) = uolr- Then, there exists a unique solution u o/O satisfying 


w GC([0,r]; H'i(D))nC^([0,r]; L2(D)), 


and there exists C > 0 such that for any t G [0, T], we have 

\\duU\\i2(Y,) + \\u{t, •)llL 2 (n) < C (||/||_H-1(S) + ll“lllL2(r2)) • 

From the above Lemma one can see that, if {uq, ui) = (0,0), the Dirichlet-to-Neumann map Aq is con¬ 
tinuous from to L^(S). Therefore we denote by ||Ag|| its norm in L^(S)). 

1.2 Main results 

In order to state our main results we first introduce some notations: 

Let r > 0 such that T > 2r and D C B{0, ^) = |x G M”, |x| < ^|. We set Qr = [0,T] x B{0, ^). 
We consider the following sets 
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Note also Q* = H Let denote by Q* = QnQ*. We remark that Q* is made of lines making 
45° with the f-axis and meeting the planes t = 0 and t = T outside Q^. We denote by Qjj = Q n 
We remark that Qj is made of lines making 45° with the f-axis and meeting only the planes t = 0 outside 
Qj.. Let’s note that Q* C Qtt C Q. 


Remark 1 In the particular case where = i?(0, ^), we remark that Q* = which is the region I in 
Figur^L^ And which is the region / U 1/ U III U IV. 



M=,_- 


'\ \^\ = T---t 


Figure 1.2 

Further, given go £ C^iQr) M > 0, we introduce 

A*{qo,M) = {qe C^iQr), g = go in Qr\Q*, Ml-°{q) < M} , 


and 

,4*'(go,M) = {g G C^{Qr), g = go in Qr\Qi, l|g||L°°(Q) < M} . 

Then our first main result can be stated as follows: 

Theorem 1 Assume that T > 2Diam (fl). Then, for every qi, q 2 G ^*(go, M), there exist two con¬ 
stants C > 0 and pi G (0,1), such that we have 

l|gi -g2||H-i(Q*) ^ (II^Qi - + |iog||Aqi - Agalir^), 

where C depends only on fl, M, T, and n. 
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Tl 

Suppose in addition that qi, q 2 £ for ■s > ^ lki||_ff‘>+i(Q) ^ * = !> 2, for 

some M > 0, then there exist two constants C" > 0 and p 2 £ (0,1) sued tdat 

Iki - d2||L“(Q.) < C*' (l|Aqi - Aqall + |log||A<ji - Agjlll"^)^". (1.2) 

As an immediate consequence of Theorem [T] we have the following uniqueness result. 

Corollary 1.1 (Uniqueness) Under the same assumptions, for every qi, q 2 G Al*((70) Af), we have the 
uniqueness 

Aqi(/) = Aq^(/), /or any/G imply qi{t,x) = q 2 {t,x), 

everywhere in Q*. 

Let us note that in this result we determine the time dependent coefficient q from full boundary mea¬ 
surements Aq only in a subset Q* C Q, provided that it is known outside of this part. 

In order to extend this result to the determination of q in a larger region D Q* we need more in¬ 
formation about the solution u. Namely we need the measures of {u{T, .),dtu(T, .)). So, let’s introduce 
the following boundary operator: 

^q : ffi(S) ^ L2(S) X H^{n) X L^{n). 

f I— {duU, u{T,.), dtu{T,.)) 

From Lemma 0 we deduce that, if (uq, ui) = (0, 0), the operator Aig is continuous from to 

L^(S) X X L^(n). We denote by \\A?q\\ its norm in LfiTf) x x . 

Then, the second result is the following: 

Theorem 2 Assume thatT > 2Diam{fl). Then, for every qi, q 2 G A^{qo,M), there exist two con¬ 
stants C > 0 and pi G (0,1), such that we have 

WQI ~ Q2\\H-^{Qff) — C -f \ — Mq^lW ), 

where C depends only on Q, M, T, and n. 

Tt 

Suppose in addition that qi, q 2 G for ■s > ^ * = 1)2, for 

some M > 0, then there exist two constants C' > 0 and p .2 £ (0,1) such that 

Ikl - d2||L-(Q„) < C {\\Mq^ -^qall + I ^Og\\^q^ • 

where C' depends on fl, M, T, and n. 

As an immediate consequence of Theorem]^ we have the following uniqueness result. 

Corollary 1.2 (Uniqueness) Under the same assumptions, for every qi, q 2 G A^{qQ, M), we have the 
uniqueness 

^qiif) = ^q 2 {f): far any f G imply qi{t,x) = q 2 {t,x), 

everywhere in Q^. 
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With zero initial data there is no hope to recover q{t, x) over the whole domain Q, even from the 
knowledge of the boundary operator Mq. However, from measurements made for all possible initial data, 
we can extend the results in Theorem[^and Theorem[^to the determination of q over the whole domain. 
We define the boundary operator 

Xq : H'i(S) X H'i(Q) X ^ x H^{n) x 

{f,uo,ui) I—> {d„u, u{T,.), dtu{T,.)) 

From Lemme [rd| we deduce that the linear operator Ig is continuous from x H^{Q) x to 

L^(S) X X L^(n). We denote by ||Xq|| its norm. 

Then, our last result can be stated as follows: 

Theorems Assume thatT > 2Diam{Q). Then, for every qi, q 2 G C^{Q), such that ||gi||Loo(Q) < M, 
for i = 1, 2. There exist two constants C > 0 and fii G (0,1), such that we have 

Ikl -92 ||h-1(Q) <C[\\Zq^ -'^q2V^ + |log||Xg, -Xq^lir^) , 
where C depends only onQ, M , T, and n. 

Suppose in addition thatqi, q 2 G for s > | and ||9i||H'>+i(Q) < M, i = 1, 2, for some 

M > 0, then there exist two constants C' > 0 and p 2 G (0,1) such that 

Ikl — < 12 \\l^{Q) < C {\\Tq^ — Xq^ || + | log ||Xg^ — Xq^\\\ . 

As an immediate consequence of Theorem|^ we have: 

Corollary 1.3 Under the same assumptions as in Theorem^ we have the uniqueness 

Xq^ =Xg 2 , imply qi{t,x) = q 2 {t,x), in Q. 

This paper is organized as follows. In section]^ we construct special optics geometrical solutions to 
the wave equation Using these geometric optics solutions, in section we prove Theorem [T] in 
section 1^ we prove Theorem]^ and in sectionwe prove Theorem]^ 

2 Geometric optics solutions 

In the present section, we collect some results which are needed in the proof of our main results. We 
start by the following Lemma (see ifT^ . pFl): 

Lemma 2.1 Let T > 0 and q G L°°{Q), suppose that F G X^(0, T; The unique solution u of 

the system 

{ {pi — A + q{t, x)) u{t, x) = F{t, x) in Q, 
u{0,x) = dtu{0,x) = 0 in U, 

u{t, x) = 0 on S, 

satisfies 

ueC{[0,T];Hpn))nC\[0,Ty,LpQ)). 

Moreover, there exists a constant C > 0 such that 

\\dtu{t, .)IIl2(q) + l|Vrt(f, 01^2(0) < C'll-^lli,i(0XX^(tt))' (2-3) 
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Using Lemma [2U] we are able to construct suitable geometrical optics solutions for our inverse problem, 
which are key ingredients to the proof of our main results. 


Let ip G Notice that for all oj G = {w G |cj| = 1}, the function 

a(f, x) = ip{x + tu) 

solves the transport equation 


{dt — uj.V)a{t, x) = 0. 


(2.4) 

(2.5) 


Let’s now prove the following Lemma: 


Lemma 2.2 Let q G C^{Q) such that < M. For oj G EF and (p G C“(M"’), we consider 


( 2 . 6 ) 


the function a defined by ( |2.4| ). Then, for A > 0, the equation 

— A + q{t, x)) u{t, x) = 0 in Q, 


admits a solution 
of the following form 
where R^{t, x) satisfies 
and 

Moreover, 


G C([0, T]; Lfi(U)) n C^lO, T]-L^{n)), 

x) = a{t, + R^{t, x), 

R^{t,x) = 0, for all {t,x) G S 

dtR'^{0,x) = R~^{0,x) = 0, X G U, 
dtR~ (T, x) = R~ (T, x) = 0, X G U. 


(2.7) 


0 , 


( 2 . 8 ) 


II^"^IIl 2(Q) + ||Vi?=^||i2(Q) < c I 

where C depends only on Q, T and M. 

Proof. We adapt the strategy developed in the proof of a similar result in ifTSl , where a time independent 
potential q was considered. In light of \2.6\ and (2.7 1 it is enough to prove the existence of R^ satisfying 


—A + q{t,x)) x)=— — A + q{t,x)) ^o(f, in Q, 

R^{e,x) = 0, dtR^{9,x) =0, 6 = 0, or T in U, (2-9) 


R^{t, x) = 0 


on S, 

and obeying ( |2.8| ). We prove the result for ri+. The existence of u~, being handled in a similar way. To 
do that note 

g{t,x) = - {d^ - A + q{t,x)) 

and use ( |2.5| ), getting 

g(^t,x) = {d^-A + q{t,x))a{t,x) = -F^^^-‘^+^'>go{t,x), (2.10) 

where go G (0, T ; (fl))- Thus, ii is a suitable solution to the system (2.91 satisfying 

RGC{[0,T]-H^mnc\[0,T]-L\n)) 
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and the function 


w{t,x) = / R{s,x)ds 

Jo 


solves the following equation 


Where 


((9| - A + q{t,x)) w{t, x) = Fi{t,x) + F 2 {t,x) in Q, 

'w{0,x) = 0, dtw{0,x) = 0 in fl, 

w{t, x) = 0 on E. 

Fi{t,x) = / g{s,x)ds, and F 2 {t,x)= / [q{t,x) — q{s,x)]R{s,x) ds. 


/o JO 

Let r G [0, T]. In use of Lemma [ aT] on the interval [0, r], there exists a constant C > 0 such that 

\\dtw{T,.)\\l,^^^ < C\\F,+F2\\ 

Ll( 0 ,T;L 2 (n)) 

- (|I^iIIl 2 (Q) + ll^ 2 |li.(o ,r;L 2 (Q))) • 


Using (|2.11|), we have 


\\^^\\h{0,r,L^n)) ^ <^T|kllioo(Q) / \\dtw{s,.)\\l2^^)ds. 


Then, it follows from p.l3| ) that 

\\dtw{T, .)lli2(Q) < c (^||i^i||i2(Q) + Ikllioo(Q) ^ || 5 tu;(s, .)lli2(o) ds 

Then, from GronwalTs inequality, one gets 

where the constant Ct > 0 depends on T and ||g'|| 2 ,oo. From where we get 

Plli2(Q)<CT||Fi||i,(Q), 

according to ( |2.11[ ). Further, as 

11-^1111-2(0) = ^5^2 / l/ gois,x)dsiF^^^-‘^"^'''>)ds\^ dxdt, 


'Q Jo 


by (2.10l and(2.12i. Then, integrating by parts with respect to s, we deduce from (2.14) that there 
a constant C > 0 such that 

(j 


Finally, Since || 5 ||l 2 (q) < C ||<y7||_H'3(RrL), using the energy estimate (2.31 for the problem (2.91 we 

II'^-^IIl2(q) < Ct ||¥’||/r3(]Rn), 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


exists 


obtain 


This completes the proof. 


□ 












3 Proof of Theorem [I] 

In the present section we will prove a log-type stability estimate in determining q appearing in the initial 
boundary value problem o with {uo,ui) = (0, 0). The main ingredients of the proof are geometric 
optics solutions introduced in Section and X-ray transform. We start by considering geometric optics 
solutions of the form (2.7 1 . We only assume that supp ip C £/r, in such a way we have 


:^n—l 


supp y? n n = 0, and (supp p±Tuj)r\Q = 9, VcuGS’' 

Then we have the following preliminary estimate which relates the differential of two potentials to the 
Dirichlet-to-Neumann map. 

Lemma 3.1 Let qi,q2 £ ■A*{qo,M), and put q = {q2 — Qi)- There exists C > 0, such that for any 
io G and p G C^{x^r)< the following estimate 


q{t,x — too) p^{x) dxdt\ < C — Ag^|| -|- 


*)• 


(3.15) 


holds true for any sufficiently large A > 0. 

Proof. In view of Lemma 2.2 and using the fact that supp (/p n = 0, there exists a geometrical optics 
solutions M 2 ,a to the equation 

{dl - A -h 92(f, x)) U 2 ,xit, x) = 0 in Q, U 2 ,x\t=o = dtU 2 ,x\t=o = 0 in If, 

of the form 


U 2 ,x{t, x) = a{t, -f R 2 ,x{t, x), 


(3.16) 


where i? 2 ,A satisfies 
and 


dtR2,Xlt=0 — R2,Xlt=0 — 0 ) R2,XIE 


= 0 . 


c 


||-R2,a||i,2(q) < — 


(3.17) 


We denote by ui, the solution of 


(df - A + qi(t,x)) ui(t,x) = 0 in Q, 
ui(0, x) = dtui(0, x) = 0 in n, 

ui(t,x) = U 2 ,x(t,x) := fx(t,x), on S. 


Putting u(t, x) = ui(t, x) — U 2 ,x(t, x), we get that 


- A -h 9 i(f, x)) u{t, x) = q{t, x)u 2 ,a(L x) in Q, 
u{0,x) = dtu{0,x) = 0 in n, 

u{t, x) = 0 on S. 


Applying Lemma 2.2 once more for A large enough and using the fact that supp p iL Tcu fl Q = 0, we 
may find a geometrical optic solution vx to the backward wave equation 

(d^ -A + qi(t, x)) vx(t, x) = 0, in Q, vx\t=T = dtVx\t=T = 0, in 
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of the form 

where Ri x satisfies 

and 


vx{t, x) = a{t, x)e + Ri^x{t, x), 

dtRi,\\t=T = Ri,\\t=T = 0, -Ri,a|s = 0, 


(3.18) 


(j 

Consequently, by integrating by parts and using the Green’s formula, we obtain 

/ q{t,x)u 2 ,x{t,x)vxit,x)dxdt = [ {d^ - A + qi{t,x))u{t,x)vxit,x)dxdt 

Jq ’ JQ 

= {Ag^ - Ag^)fxit,x)vxit, x) da dt, 


(3.19) 


(3.20) 


So, ( |3.16| ), ( |3.18| ) and ( |3.20| ) yield 


/ q{t,x) a‘^{t,x) dx dt + / q{t,x)Ri^x{t,x)R 2 ,xit,x) dx dt 
IQ Jq ’ ’ 


+ 


f q{t,x)a{t,x) dxdt 

J Q 


'Q 

(Aga - Agi)/A(f,3:) Vx{t,x) da dt. 


(3.21) 


From ( |3.21| ), ( |3.17[ ) and ( |3.19| ) it follows that 


\ j^q{t,x)o?{t,x)dxdt\ < {{Ag^ - Ag^)fx{t,x)vx{t,x)\ da dt + j ||(^||^3(]Rn), 

where the constant C > 0 does not depend on A. Hence from the Cauchy-Schwartz inequality and using 
the fact that fx{t, x) = U 2 ,x{t, x) on S, we obtain 


f C 

/ q{t,x)a^{t,x)dxdt\ < ||A ,2 - A^iH \\u2,\\\hRi:) II^^a||l 2 (s) + — 
J Q 


O’ 


(3.22) 


Further, as Ri^x\T, = 0> for f = 1, 2, we deduce from ( 3.22[ ) that 

q{t,x)a?{t,x)dxdt\ < C |^||Ag 2 - AgJ| ||m 2 ,a - R 2 Ah^{q) II^^a - i(i,A||/ri(Q) + ^ 


'Q 


Bearing in mind that 


I^A - -Ri,A||jri(Q) < C\ ||(/5||/r3(]Rn), 


we end up getting that 


\\U2,X — R2,\\\h^{Q) < CX^ 

J q{t,x) 0 ?{t,x) dxdt\ < C - AgJ| + ^ 


O' 


Therefore by extending q{x, t) by zero outside Qr and recalling (2.4 1 , we find out that 


q{t, X — tui) (p^{x) dx dt\ < C [ A^||Aq 2 — Ag^ || + ^ 


A 


This completes the proof of the Lemma. 


□ 
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3.1 X-ray transform 

The X-ray transform R maps a function in into the set of its line integrals. More precisely, if 

u G and {t, x) G 

:= [ f{t,x -tu})dt, 

Jr 

is the integral of / over the lines {{t,x — tui), f G M}. 


Using the above Lemma, we can estimate the X-ray transform of the differential of potentials as fol¬ 
lows: 

Lemma 3.2 There exists a constant C > 0, /3 > 0, 5 > 0, and Aq > 0 such that for all oj G we 
have 


\R{q){uj,y)\<C 


'■Ql 


1 


a.e. y G 


for any A > Aq. 

Proof. Let f gCo“(M") be a positive function which is supported in the unit ball 5(0,1) such that 

II<?^IIl 2(R") = 1- Define 

iPe{x) = 

where y £ £/r- Then for sufficiently small e > 0 we can verify that 

supp ipe nTl = $, and supp ip^ ± Toj n U = 0. 


And we have 

rT 


/ q {t,y — toj) dt\ = \ / / q{t,y — tu}) (p‘^{x) dx dt\ 

Jo ^ Jo Jr^ ^ 


< 


/ / q{t,x — tuj) ipl{x) dx dt\ + \ / / {q{t,y — too) — q{t,x — tcj)) (p‘^{x) dx dt\. 

Jo dK" Jo Jr^ 

Since Ugllc^Q) — we have 

\q{t, y — too) — q{t, x — tLo)\ < C \x — y\. 

Applying Lemma [3T] with p = we obtain 

1^ q{t,y-tuj)dt\ < C (^A^IIAqa - AgJI -f ^ + C \x - y\p‘^^{x) dx. (3.23) 

On the other hand, we have 


¥’£||H3(Rn) < Ce / \x-y\pi{x)dx <Ce. 


Thus, from ( |3.23[ ), we obtain 

rT 


'0 


q{t,y - tu})dt\ < C |^A^||A ,2 - AgJ| -f e ® -h Ce. 
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We select e such that 


--e 


e = 

Then there exist constants 5 > 0 and /3 > 0 such that 

rT 


X 


q{t, y - tui) dt\ < C || - AqJ| + 

Using the fact that q = 0, outside Qr, we get 

I J q{t,y - toj) dt\ < C ^A^||Aq 2 - AqJ| + , a.e. y G w G S' 


n—1 


On the other hand, if |y| < then 


Indeed, we have 


y(f ,y — toj) = Oyt G 


\y - tuj\ > |f| - |y| > f - 


(3.24) 


(3.25) 


(3.26) 


So that, if f > -, from (3.261, we have {t, y — tuj) ^ . And if f < -, we have also (f, y — tuj) ^ . 

Consequently, 


(f, y — tuj) ^ D Q*, for all f G M. 


Using the fact that q = q 2 — qi = ^ outside Q*, we deduce (3.25 1 . Therefore, 

[ q{t,y-tui)dt = 0, a.e. y G 5(0, ^). 

Jr ^ 

By a similar way, we prove that in the case where |y| > T — we have 

(f, y — tuj) ^ D (3*, for all f G M. 

Then we conclude that 

/ q{t,y — toj) dt = 0, a.e. y ^ jz/r, tu G S"'“^. 

Jr 

Consequently, by ( |3.24 i and (3.27 1 , one gets 

l^(7)(w,y)| = \ j^q{t,y -tu})dt\ < C (^A^IjAg^ - AgJ| + , a.e. y 

This completes the proof of the Lemma. 

Let now 

5 = {(r,OeMxM", |r|<|e|}, 

and let the Fourier transform of q G L^{R^+^) 


(3.27) 


GM'^, wGS’"-^ 


□ 


q{r,0 = 

Our goal now is to prove the following 


q{x,t)e ^^'^dxdt. 
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Lemma 3.3 There exist constants C'>0, /3>0,5>0 and Aq > 0 such that the following estimate 
holds 


T,0\<C 


91 


1 


for any (r, ^ E and A > Aq- 

Proof . Let {t,^) £ E and ( G 


such that = 0. By defining 




we have a; G S"" ^ and uif = r. 


By the change of variable x = y — tco we have for all ^ G w G ^ 

[ R{q){uj,y)e~"y-^ dy = [ ( f q{t,y - tu) dt^ dy 

Jr J 


= 

Denote (r, = (uifj) G E. Since supp q{t, .) C D C -6(0, £), then we have 


/ 

JR^ 


"nB(0,§+T) 

In terms of Lemma the proof is completed. 


R{q){^,y)e dy = q{Tj). 


□ 


3.2 Stability estimate 

We are now in position to complete the proof of Theorem[^ For p > 0 and 7 G (NU {O})”"''^, we denote 
I 7 I =71 + ... + 7n+i, 6(0, p) = {x G |x| < p}. 

We consider the following Lemma 

Lemma 3.4 (see 4251/ ) Let O be an open set o/6(0,1), and F an analytic function in 6(0,2), satisfying 
the following property: there exist constant M, p > 0 such that 

I|8"^’IIi~(B(0.2)) < V7 £ (N U {0})”+'. 

Then, 

II-^IIl°°(s(o,i)) < (2M)^ ^II-^IIl°°(0)- 

where p G (0,1) depends on n, p and \0\. 
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The Lemma is conditional stability for the analytic continuation, and see Lavrent’ev, Romanov and 
Shishat-sKii. ifTdll for classical results. For fixed a > 0, let us set 

= q{a{T, 0 ) for (r,^ € 

It is easily seen that Fa is analytic and we have 

\d'^Fa{T,i)\ = \d'^q{a{T,i))\ = \d'^ [ 

7r"+i 

(3.28) 


jn + l 


Therefore, from (3.281 one gets 

|5^Fa(r,0l < [ |g(f,x)|al'^l(|x|2 + f2)^dxdf < ||g|Ui(Q^) {2T^)'-^ < C 

7r"+i (2 


)l7l 


e . 


Then, applying Lemma 3.4 in the set (9 = S H B{0, 1) with M = (9e“ and rj = T we can take a 


constant fj, € (0,1) such that 

|i"a(r,OI = |9 («(t, 0)I < (r,0 G 5(0,1). 

Hence, by using the fact that aE = {a(r, ^), (r, £ E} = E, we get for (r, |) G 5(0, a) 


(r,e)| = |5„(«-i(r,0l < ||5„||(( 


L°°(0) 




L°°(S( 0 ,a)nE) 


L°°{Ey 


(3.29) 


On the other hand we have 


|2//i 

lj^-l(Rn + l) 


= / (l + l(r,0Pr'l9(r,0p7rde + / {I + \{r,0\^)-^\q{r,0\^ drdA 

\7|{t,^)|<o 7|(r,OI>a / 

< c(a^+^ ll9lli-(S(0,a)) + ll2lli2(R«+l)) 


\ i/ri 


From (3.291 and applying Lemma 3.3 we obtain 


/ 1 \ 

< C U"+^e2“(i-'^)(A^||A,, - AgJI + ^)2M + a-2j 

/ n+l 2a(l-<j) „„ „ rt + 1 2^(1-^) „c r,/,, \ 

< C I a e A^^||Aq 2 — Ag^|p + a e A + a /^ | . 


|2//i 

Ih-i(]R"+A 


Let ao > 0 be sufficiently large and a > aq. Set 


n+3 oc{l — iJ.) 

A = a 2/^15 e . 


By a > ao, we can assume that A > Aq, and we have 


n + l 2a(l /i) ^ j. ^ / 

A“2^ = a“2/^. 
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Then 


I 2/m 


(5(n+l)+;3(n+3) 2a{5+^){l —^i) 

< C i a e 


< 


11^92 - V IP + 

C ^e^"||Ag2 - Agjp + , 


where N depends on 6, /3, n, and [i. In order to minimize the right hand-side with respect to a, we set 

1 


where we assume that 


« = ^1 l0g||Ag2 -Aq,\\ I, 


0 < II 


It follows that 


kl|jf-l(Q*) < Ikll_ft-1(R’^+1) < C* (11^2 “ ^91II + I log 11^2 “ ^91II I 

< C(||A,,-A,,r/2 + |log||A,,-A,J||-i). 


m/2 


The estimate (1.2), is now an easy consequence of the Sobolev embedding theorem and an interpolation 
inequality. Let 5' > 0 such that s = n/2 + 26'. Then, we have 


\l°°{q+ < C\\q\\Hs(Q+ 

Ih-P(Q.) II^II^'+MQ*) 


< c 


for some (3 G (0,1). Then the proof of Theorem[^is completed. 


4 Proof of Theorem |2] 


This section is devoted to the proof of Theorem We will extend the stability estimate ( |1.2| ) given in 
TheoremjTJ to an estimate in a larger region D Q*. Differently to Theoremj^ here the observations 
are given by the boundary operator introduced in Subsection |1.2[ We need to consider geometric 
optics solutions similar to the one used in the previous section, but this time, we will only assume that 
supp y? n D = 0. (We don’t need to assume that supp + ± Tui n D = 0). Let’s first recall the definition 
of the operator Mq'. 

^q : ^ L2(S) X X L2(Q). 

/ I—> {d^u, u{T,.), dtu{T,.)). 


We denote by 


if) = dyUj, if) = UjfT ,.), (/) = dtUj{T ,.), for y = 1, 2. 


Lemma 4.1 Let qi,q2 G -A^iQo, M), ip G such that supp <y9 n D = 0, and put q = iq2 — <?i)- 

Then, there exists C > 0, such that for any uj G the following estimate 

f q{t,x-toj)p‘^{x)dxdt\ <C (^X^\\Mq 2 -fgqf\ + ^ (4.30) 

holds true for any sufficiently large A > 0. 
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Proof. In view of Lemma 2.2 and using the fact that supp (/p n 12 = 0, there exists a geometrical optics 
solutions M 2 ,a to the equation 


{dl - A + 52(f, x)) U 2 ,xit, x) = 0 in Q, U 2 ,x\t=o = dtU 2 ,x\t=o = 0 in 12, 


of the form 


where i? 2 ,A satisfies 
and 


U 2 ,x{t,x) = +R 2 ,x{t,x), 


dtR2,X\t=0 — R2,X\t=0 — 0) -R2,A|S — 0, 


(4.31) 


C 


I-^2,a||i,2(q) < — ||v^||/r3(]Rn) 


(4.32) 


We denote by ui^x-> the solution of 


{df - A + qi{t,x))ui^x{t,x) = 0 in Q, 
ui,x{0, x) = dtui^xi^, x) = 0 in 12, 

x) = U 2 ,x{t, x) := fx{t, x), on S. 

Putting ux{t, x) = x) — U 2 ^x{t, x), we get that 

(Of - A + 5 i(f, x)) ux{t,x) = q{t,x)u 2 ,x{t,x) in Q 
ua(0,x) = dtux{0,x) = 0 in 12 

ux{t, x) = 0 on S. 

Applying Lemma [2Al once more for A large enough, we may find a geometrical optic solution vx to the 
backward wave equation 

(of - A + 5 i(f,x)) XA(f,x) = 0, in Q, 


of the form 

where Ri x satisfies 
and 


vxit, x) = a{t, x)e + Ri^x{t, x), 


dtRi^x\t=T — Ri,x\t=T — 0, -Ri,a|s — 0, 


C 


(4.33) 


11-^1,a||l 2(Q) < — ||(p||fr3(]Rn). 

Consequently, by integrating by parts and using the Green’s formula we obtain 


(4.34) 


'Q 


q{t, x)u 2 ,xit, x)vx{t,x) dx dt = 


)(fx)vx(t,x)dadt 


"92 ‘^'"91 


+ / {Ml-^l){fx)dtvx{T,.)dx 
Jn 

- f «-^3j(/A)xA(r,.)dx, 

Jn 


(4.35) 
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Then, by replacing U 2 ^\ and v\ by their expressions in the left hand side of (4.351 and using (4.321 and 


(4.34|), then from Cauchy-Schwartz inequality, one gets the following estimate 


C 


'Q 


q{t,x)a‘^{t,x)dxdt\ < ||(^g\ - ^g\)(/A)||L2(s)lkA||i,2(s) + y 

+ ll«-0(/A)llL2(n)l|5iUA(r,.)||i.(n) 

+ ll«-O(/A)llL2(0)lbA(T,.)||i2(^). 


Then we obtain. 


f q{t,x)a^{t,x)dxdt\ < (||(^g 2 -^g\)(/A)|li 2 (E) + \\i^q 2 - ^qi)ifx)\\miQ) + \\i^q 2 - 
Q 


Setting 


^A|li2(s) + lkA( 2 ', •)IIl2(q) + \\dtvx{T, 


</>A = (t'A|E> vx{T ,.), dtvx{T ,.)) 


C 

+ -r 


(4.36) 


Then, from (4.361, we get 


C 


'Q 


q{t,x)a {t,x)dxdt\ < \\{^q 2 —^qi){f\)\\L^(T)xm(n)xL^{n) ll'/’A||L 2 (s)xL 2 (Q)xL 2 (r 2 ) + y||‘f'il/r3(]Rn) 


and using the fact that fx{t, x) = U 2 ^x{'t, x) on S, we obtain 


C 


q(t,x)a (t,x)dxdtl < - ^q^W |k 2 ,A|bl(S) ll</'A||L 2 (s)xL 2 (n)xL 2 (n) + Y 


JQ 

Further, as Ri^x\s = 0> for * = 1)2, we deduce that 

\ j^q{t,x)a‘^{t,x)dxdt\ < C - ^qiW ||** 2 ,A - R2,x\\h^q) Ui,x\\m{Q)xLHn)xL^in) + j 

where 

</'i,A = (**A - Ri,x, vx{T ,.), dtVx{T ,.)). 

Using the fact that Ri^x{T, •) = dtRi,x{T ,.) = 0 on U, we have 

\\u2 - R2 \\h'2{q) < ll¥’ll/r3(]Rn), 

and 

II</’1,a||h1(Q)xL2(q)xL2(Q) < I^A - Ri,x\\h^{Q) + l|r'A|t=TllL2(fj) + \\dtVx\t=T\\L^{n) 


Therefore by extending q{t, x) by zero outside Qr and recalling (2.4 1 , we find out that 


I / / q{t^x — tuj) q^‘^{x) dx dt\ < C [ X 

Jo Jr" 

This completes the proof of the Lemma. 

Let’s move now to prove the following Lemma 


cq2 ^oq^ 


+ - 


‘)- 


□ 
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Lemma 4.2 There exists a constant C > 0, /3 > 0, 5 > 0, and Aq > 0 such that for all w G S 
have 

\R{q){uj,y)\ < C -^qiW + , a-e- y G K”- 

for any A > Aq. 


we 


Proof . We consider {ipe)e defined in the proof of Lemma 3.2 We only assume that y ^ Tl, then for 
sufficiently small e > 0, we can verify that supp n = 0. Taking in acount this last remark, using 
Lemma 4.1 and repeating the arguments used in Lemma [3]^ we obtain this estimate 


\ j^q{t,y -tui)dt\ < C (^X^\\^q 2 -^qiW + > a.e.y ^ B{0, ^). 

On the other hand, if y G B{0, ^), then we have 

q{t,y — tu) = 0, VtGM. 

Indeed, we have 

\y-tuj\ > |f| - |y| > t - 


(4.37) 


(4.38) 


(4.39) 


So that, from (4.391, we deduce that for all f > | we have (t, y — tui) ^ And if f < |, we have also 
that {t, y — tui) ^ We recall that = Q H TXf. Consequently, we have 

(t, y — tuj) fi Q((, for all f G M. 


q{t,y -tuj)dt = 0, a.e. y G 5(0,-). 


Then, using the fact that y = ^2 — yi = 0 outside Qjj, we obtain (4.381. Therefore 

(4.40) 

z 

In light of (4.37 1 and (4.401, the proof of Lemma [4^ is completed. □ 

Using the above result and in the same way as in Sectionj^ we complete the proof of Theorem]^ 


5 Proof of Theorem |3] 

In this section we deal with the same problem treated in Section and except our data will be the 
response of the medium for all possible initial data. As usual, we will prove Theorem [fusing geometric 
optics solutions constructed in Section and X-ray transform. Let’s first recall the definition of the 
operator Xq! 

Iq : 5'1(S) X H^{n) X X2(U) ^ L2(S) x H^{n) x L‘^{n). 

V' = if, Uo, Ui) I—> {d,yU, u{T ,.), dtu{T ,.)). 

We denote by 

= duUj, = Uj{T, .), Xj(V’) = dtUj{T, .), for j = 1, 2. 
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Lemma 5.1 Let qi, q 2 ^ C.^{Q)> cmd put q = {q 2 — qi)- There exists C'>0, /3>0, 5>0 and Aq > 0 
such that for any uj G we have the following estimate 


\R{q){uj,y)\<C 




a.e. y G M"'. 


for any A > Aq- 


Proof. Let p G For A sufficiently large, Lemma 2.2 guarantees the existence of the geometri¬ 

cal optics solution U 2 ^\ to 


of the form 

where i? 2 ,A satisfies 

and 


{d^ -/\ +q 2 {t,x))u 2 ,\{t,x) = Q, in Q, 


x) = a{t, -f i? 2 ,A(f, x) 


9tTl2,X\t=0 — -R2,A|f=0 — 0) -R2,A|S — 0) 


(5.41) 


c. 


1^2,\\\l2{Q) < Yll‘^ll/t3(]Rn) 


(5.42) 


We denofe ui^\ fhe solufion of 


{df - A + qi{t,x))ui^x{t,x) = 0 in Q, 

rii,A(0, x) = U 2 ,x{ 0 , x), dtUi^xiO, x) = 9*^2,a(0, x) in Q, 
ui,x{t,x) = U 2 ,x{t,x) := fx{t,x), on S. 

Puffing ux{t, x) = x) — M2 ,a(L x), we gef fhaf 

- A + qi(t,x)) ux(t,x) = q(t,x)u 2 ,x(t,x) in Q 
ux(0,x) = dtux(0,x) = 0 in n 

ux(t, x) = 0 on S. 

Applying Lemma once more for A large enough, we may find a geomefrical optic solufion vx to the 
backward wave equation 

(di - A +qi(t,x)) vx(t,x) = 0, in Q, 


of the form 

where Ri x satisfies 

and 


vx(t, x) = a(t, x)e + Ri^x{t, x), 


dtRi,x\t=T — Tti^x\t=T — 0) -^i,A|s — 0, 


C 


(5.43) 


11-^1,Alli2(Q) < — ||(^||fr3(]Rn). 

By infegrafing by parfs and using fhe Green’s formula, one gefs 


(5.44) 


IQ 


qf, x) U 2 ,x{t, x)vxit, x) dx dt = j - Zl^){fx)vx{tix) da dt 


19 



+ [ {Xl - II) (V^a) dtvx iT,.)dx- [ {Zl - ) (V^a) (T, .)dx, (5.45) 

Jo. Jfl 

where 

V’A = ('W2,A|E) ^2,A|t=0) ^t^^2,A|t=o)- 

Next, we proceed by a similar way as in the proof of Lemma [4~T| we get 


I / q{t,x)a^{t,x)dxdt\ < WZq^ - Zq^\\ \\'ij;x\\mCE)xH'^{n)xL'^(n)\\4>x\\L'^(j:)xL'^{n)xL^(n) 
•JQ ^ 

where 

</>A = {vx\^, Vx\t=T, dtVx\t=T)- 
Further, as i?j_A|s = for * = 1) 2„ we deduce that 


q{t,x) a‘^{t,x) dx dt\ < 


11^92 “^91II ll^l,A||H2(Q)x/tl(r2)xL2(Q)ll<?^l,A||_ffl(Q)xL2(n)xL2(n) 

+ y ll‘^llH3(Rn), 


where 

4>1,\ = {vx - Rl,\, Vx\t=T, dtVx\t=T) , 'lpl,X = {U 2 ,X “ f? 2 ,A, «2,A|t=0, dtU2,x\t=o) • 
Using the fact that Ri^x{T ,.) = dtRi^xiZ ,.) = 0 on U, we have 

ll<?^l,A||_ffl(Q)xL2(f7)xL2(Q) < CA ||(/?||jy3(]Rn), 

On the other hand, since i?2,A(0, ■) = dtR2,x{0 ,.) = 0 on Q, we have 


^l,A||_H'2(Q)xHi(ft)xL2(n) < 1^2,A — f? 2 ,A||_ff 2 (Q) + ||rt2,A|t=olli^i(n) + l|5trf2,A|t=ollL2(f2) 


Therefore by extending q{t, x) by zero outside Q and recalling (2.4 1 , we find out that 

I f [ q{t^x — toj) qx^{x) dx dt\ < C ()?\\Zq^—Zq^W 
Jo Jr" \ 

Now, in order to complete the proof of Lemma 5.1 it will be enough to fix y G consider {ipe 


defined as before, and proceed as in fhe proof of Lemma [T2| By repeafing fhe argumenfs used in fhe 
previous secfions, we complefe fhe proof of Theorem]^ □ 
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